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ABSTRACT 


Individual massive stars with M^ ^ <-6 have huge stellar winds that 
create interstellar bubbles. Stars wx?h masses greater than 8M© (4M0?) 
are considered supernova progenitors. These massive stars are numerous in 
OB associations where few supernova remnants are detected. Model calculations 
describing the evolution of an association show: 1) that layge, hot cavities 

are formed by pushing the ambient gas into neutral shells; 2) that the shell 
radii change with galactocentric radius, 3) that only thirty percent of the 
interstellar medium is in the form of supercavities and 4) that a consequence 
is that only a small fraction of supernovae form supernova remnants. These 
results have strong bearing on interpretation of interstellar studies being 
done by lUE and by HEAO-B. 

INTRODUCTION 

i 

The previous talk (ref. 1) described the observable interstellar 
bubble being driven by an 0 star, HD148937. Other talks (ref. 2,3) discussed 
lUE observations of supernova remnants which presumably had massive stars as 
progenitors. We know that O and B stars form as associations throughout 
the galactic plane. If individual stars can dramatically alter the inter- 
stellar medixim, we are led to wonder what an OB association would do. The 
model calculations discussed here described the bulk properties of how the 
surrounding interstellar medium is changed as an OB association evolved. The 
resultant structure of the interstellar mediimi ^an be studied by lUE and 
HEAO-B and indeed such studies are underway at present. 

THE BASIS FOR CONSIDERING SUPERBUBBLES 


Our interest in the structure of the interstellar medium has been 
whetted by a recent emission line survey of the Milky Way (ref. 4) . It 
presents a photographic record of nearby dust clouds and ionized hydrogen 
regions with a rough excitation classification being possible by singly- 
ionized sulfur imagery and by doubly- ionized oxygen imagery. Several new 
supernova remnants (ref. 5, 6) were discovered from this data. Many inter' 
stellar bubbles in the form of bowshocks (ref. 7), arcs and shells (ref. 8 
and 9) were newly detected. However, ionization structure around OB 
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associations varied considerably. We were led to ask how a supernova 
remnant, formed within an OB association, would be detected. The answer 
proved non-trivial, but opened up the realization that OB associations play 
a major role in determining structure of the interstellar medium. 


Supernova remnants within our galaxy present an enigma; there do not 
seem to be enough. Based upon statistics of supernova events within similar 
galaxies, we should expect one every thirty years (ref. 10). If we assume 
that pulsars are stellar remnants left over from supernova events, then pulsiinr 
statistics imply a supernova every ten years (ref. 11). We might expect that 
for every supernova there should exist, for a finite time, a supernova 
remnant generated by the ejecta propagating through the interstellar medium. 
Yet the supernova rate derived from statistics of supernova remnants is 
about one every eighty years (ref. 12 ). 

Related to this puzzle is the apparent lack of supernova remnants 
within OB associations. There simply is a lack of non-thermal radio emission 
plus associated filamentary optical shells within OB associations. 

An Intriguing clue comes from 21-cm studies (ref. 13), that large shells 
of neutral hydrogen having 100 to 2000 pc radii were scattered throughout the 
galactic plane. Because the associated kinetic energies exceeded the 
available energy of a supernova by a factor of 100, Heiles (ref. 13) 
postulated a type III supernova event. We find a more reasonable explanation 
is available based upon evolutionary products of an OB association. 

AN EVOLUTIONARY MODEL OF OB ASSOCIATIONS 


A simple model was devised based upon the kinematic properties of stellar 
winds plus supernovae in an OB association (ref. 14). The model could be 
divided into three phases; A. Bubble Thase, B. First Supernova Phase, and 
C. Late Supernova Phase. The bubble phase lasts for about 3x10^ years during 
which the approximately thirty stars with <-6 push the interstellar 

medium away with the combined stellar winds. For a typical OB ass£ciation 
these most massive stars have the following averaged properties; = 78 . 8 , 

M~10-6Mg, yr“l, v^.~2000 km s"!, where Mbd is the bolometric magnitude, M 
is the mass loss rate in solar masses, M^, per year, and v^ is the terminal 
wind velocity. For an individual star an interstellar bubble is driven 
outward by the stellar wind. The bubble radius is (ref. 15); 


R - 97 T 1/5 . 3/5 

R^ - 27 n L 3 g t^ pc 

with outward moving velocity of 

^ - 1/5 1/5 ^ - 2/5 . -1 

ih n Li^ ^ Km s 


where 


T 1 • ,-36 -1 

= -2 M v^ 10 ergs s 


(n = ambient interstellar medium number density and tg 


( 1 ) 

( 2 ) 

(3) 

time in lO^years) . 
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The first supernova phase begins at about 3x10 years when the 30 most 
massive stars (M>15%) become supernovae. The ejecta initially has little 
interstellar gas to shock until it has expanded to the shell radius pc) . 

Because the shell is mostly HI and H 2 gas and because the ejecta is very 
dilute when it shocks the shell, the temperature is not very high and the 
shell evolution is well described by the snowplow model for supernovae. The 
shell expands with time as: 


Rg(t) = 


“l\ 1/4 

+ 3 ^ ' 


(4) 


with 


3M V 

^ 4Mm^nir 


This phase lasts until the shell stalls or until the second wave of 
supernovae occur. 


The second supernova phase commences as stars down to 8M^ become supernovae 
This occurs at t~107 years and serves to maintain and/or enlarge the shell. 
Approximately 180 stars exceed and contribute to the second phase. The 
shell radius becomes 170 to 700 pc before stalling. We note that stars 
as low as 4 Mq become supernovae, then the shell will continue to expand 
and be maintained longer. 


The model calculations included number densities with scale heights 
and gravitational restoring forces depending upon the z— distance and for 
three galactocentric radii: 5 kpc, 10 kpc and 20 kpc. As the gravitational 

force only affects the bubble in the z-direction, bubble radii, etc. were 
calculated for parallel to the plane and perpendicular to the plane 
(Table 1). The radii perpendicular to the plane, being least affected by 
galactic rotation, are felt more reliable. We find that the bubble radius 
changes with galactocentric radius as calculated by our model (Figure 1). 

Of about 50 hydrogen shells, we find only six to be assoclable with OB 
associations. This is not surprising as the larger bubbles would surround 
very evolved associations with remaining spectral types trending towards A. 
Moreover, interstellar extinction prevents detecting the more distant 
associations. Approximately thirty percent, not ninety percent, of the 
interstellar medium is superbubble interiors. This indeed will produce only 
arcs, or partial shells as shells break into adjacent shells, setting up hot, 
ioniz ed tunnel s . 


MODELLING THE SUPERNOVA REMNANT OCCURRENCES 


We realized that supernovae occurring within an interstellar bubble 
or superbubble would not create the classically detectable supernova remnants. 
Basically, a supernova remnant is the product of ejecta interacting with the 
ambient interstellar gas (ref. 16). Until the supernova ejecta (~ Sl^j) has 
swept an equivalent amount of material, the ejecta is in free expansion. 
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By the time the supernova ejecta has expanded to encounter the shell boundary 
it is too diffuse and only a brief flash (like Barnard’s Loo^)) is detected. 
Studies, of detected supernova remnants both in the Milky Way and the 
Magellanic Clouds suggest supernova remnants are detectable only if the 
ambient density exceeds 0.1 cm"3. 

A compensating factor Includes the knowledge that many runaway stars 
occur and some B stars diffuse from the association beyond the hot cavity 
Into the neutral shell. About two-thirds of all stars are In binary systems. 
In a massive binary system, the more massive star explodes and tends to flitvj 
the less massive star away. The lower— massed star then becomes a runaway 
star. (Twenty percent of all OB stars appear to be runaways (ref. 17)), 

Those runaway OB stars that escape the hot cavity create the detectable 
supernova remnants. 

Only a small fraction escape the supercavities. Since the shell radii 
are much smaller at 5 kpc, more escape the hot cavity than at 10 kpc. No 
OB stars escape the hot cavity at 20 kpc. Hence, the model demonstrates a 
very galactoqentrlc distribution of supernova remnants. We also note that 
the OB stars that escape are on the lower mass portion of the OB star mass 
function. There Indeed may be a narrow mass range of OB stars that produce 
detectable supernova remnants. Based upon thirty percent of the Interstellar 
medium being supercavltles, the following percentages of supernovae create 
detectable supernova remnants; 

R gal = 5 kpc 15 percent to 30 percent 

= 10 kpc 9 percent to 23 percent 

* 20 kpc 0 percent 

The distribution of supernova remnants is maximum near the galactic center and 
zero beyond 20 parsecs. Radio' observations bear' this out. 

CLOSING COMMENTS 

The superbubbles, along with hot supercavities around OB associations 
provide a model for comparing Interstellar column densities, velocity 
components and degree of ionization. Perhaps the most Important test region*; 
3X& young OB associations, like Orion OB 1, where the first supernova events 
should have occurred. Very highly ionized gases should be interior, with 
rapidly expanding HII, HI and H 2 shell components. 

We note that a Crab-like supernova remnant, namely ejecta still 
Interacting with the expelled magnetic fields, will be detectable for several 
thousand years. With a 50,000 year interval between supernovae within a 
specific OB association, about six percent of the OB associations would have 
a detectable supernova remnant at any given time. However, very hot interiors 
with Isothermal X-ray emission should occur for much longer times. We 
point out that the Carina Nebula has a diffuse component seen by HEAO-B 
(ref. 18) and suggest more detailed studies be done to check our models. We 
are already Investigating several superbubbles in our galaxy using lUE and 
HEAO-B, and we are studying shell structures in the Magellanic Clouds using 
HEAO-B. 
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TABLE . EXPANSION OF SHELLS AT SELECTED DISTANCES FROM THE 
GALACTIC NUCLEUS, Rg^j^, IN THE GALACTIC PUNE 


p (kpc) 
"gal 

v(kro s*') 

Parallel 

to the Plane 

Perpendicular to the Plane 

t(ioV)* 

Rs(pc) 

n(cm*^) 

t(ioV)* 

Rjtpc) 

n(cm’^)5 



A. End of Bubble Phase 





5 

17 

3 

85 

3 

3 

85 

3 

10 

21 

3 

106 

1 

3 

105 

1 

20 

33 

3 

168 

0.1 

3 

168 

0.1 



B. End of First Supernovh Burst Phase (H fclSM^ 



5 

5 

8a6 

137 

3 

9 

170* 

0.44 

10 

5 

11 

185 

1 

10 

207 

0.37 

20 

5 

19 

357 

0.1 

17 

384 

0.054 



C. End of Second Supernova Burst Phase (H a 8M, 



5 

- 

14 

179 


- 

o 

0.44 

10 

5 

19 

251 

1 

37 

497 

0.086 

20 

* ^ 4a. ♦kail Ibakd 

5 

43 

520 

0.1 

46 

693 

0.033 


♦ Expansion has exceeded R^^i^ end the shell size Is now limited by gravitational deceleration (see text). 

{ n Is the average of the region between the Initial z and final z for each phase and Is derived as: 
ttn^ |.*2 z^ e'*/** dz 

R s : — 

4/3..(z3 - z 3 1 
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• STATIONARY SHELLS 

A EXPANDING SHELLS a? 

• MODEL VALUES 



Rgai (kpc) 


Figure 1 
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